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FOREWORD 


This volume contains the User Documentation for the NASTRAN 
3-D Analysis and the companion MESHGEN Data Generator Program. 
Each major section contains the dascriptions of the input data as 
well as many examples of using the programs. 
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6.0 NASTRAN 3-D HYDROELASTIC ANALYSIS USER MANUAL 


6.1 INTRODUCTION 

The three-dimensional hydroelasticity capability installed in NASTRAN 
allows the solution of problems involving interacting, arbitrarily-shaped 
structures and fluids. The program is intended for the vibration analysis 
of fluid-filled tanks in an acceleration field where the fluid motions 
interact with the structure displacements. Both free surface sloshing modes 
and higher frequency coupled modes may be obtained from the analysis. 

The method used to formulate the fluid/structure equations is described 
in the updates to the NASTRAN Theoretical Manual. The basis for defining the 
fluid are three-dimensional finite elements connected to fluid grid points 
defining the Eulerian pressure at a point fixed in space. The use of a 
pressure single degree of freedom at each point rather than three displace- 
ments allows a finer mesh of elements with a reasonable matrix order. 

In the formulation of the fluid/structure system the interior fluid 
degrees of freedom are transformed and removed from the solution matrices. 

The eigenvalues of the combination are extracted from small, fully dense, 
symmetric mass and stiffness matrices, efficiently processed with the "Givens" 
method. The solution matrices are defined only by the free surface displace- 
ments and the reduced structure coordinates. 

All NASTRAN modeling options are available for the definition of the 
structure. All existing options for NASTRAN Executive Control and Case 
Control data for normal modes analysis are also available for the hydro- 
elastic problems. In addition to the normal NASTRAN data, a hydroelastic 
problem requires the addition of a finite element fluid model, the specifica- 
tion of its boundaries, and the addition of special control data. The set-up 
of the hydroelastic NASTRAN data decks is illustrated in the examples at the 
end of this document. 
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For hydroelastic analysis the fluid is modeled with three-dimensional 
finite elements having shapes defined by tetrahedra (FTETRA) , wedge (FWEDGE) 
and hexagonal (FHEXl or FHEX2) volumes. Tlie fluid is assumed to be locally 
incompressible and non-visious with small motions relative to the overall free 
body displacements of the system. The following options are provided for 
defining the fluid boundary conditions. 

1. The default boundary is a rigid wall. 

2. Pure free surfaces are defined with single point constraints. 

3. Free surfaces with gravity effects are specified with CFFREE data 
cards. 

4. Fluid/structure boundaries are defined by CFLSTR data cards. 

Several alternate paths are available for the execution of the problem 

and formulation of the solution equations. These are: 

1. Direct versus Modal Structure Formulation 

In the "direct" formulation the solution matrices are defined by the 
structure degrees of freedom (after constrained and omitted points are 
removed) plus one degree of freedom for each free surface point defined on 
CFFREE data. The alternate "modal" formulation calculates the modes of the 
empty structure and uses the generalized displacements of these modes with 
the free surface degrees of freedom in the solution matrix formulation. 
Although the modal formulation requires the additional cost of another eigen- 
value extraction process, the combination system matrices will be smaller. 

This method is recommended for problems where several different fluid models 
are used with the same structure model. The structure modes need only be 
calculated once. Different fluid models may be analyzed using the NASTRAN 
restart procedure to recover the structure mode data. 
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2. Compressibility Options 


Two methods are provided for defining the compressible fluid effects. 

The overall compressibility of the enclosed volume may be specified as a 
parametric number which, in effect, provides a stiffness factor applied to 
the total volume change. The alternate method produces zero volume change 
by automatically constraining one degree of freedom in the system. The 
latter method ib not allowed in the model formulation option. 

3. Differential Stiffness Effects (Ullage Pressure) 

An option has been provided for Including the effects of ullage pressure 
on the structure stiffness. These additional stiffness terms are calculated 
in a separate structure-only Rigid Format 4 analy -is with pressures defined 
by static loads. The differential stiffness is transferred to the hydro- 
elastic problem with the NASTRAN checkpoint/restart procedure and is con- 
trolled by two parameters, DISTIF and DIFSCALE. 

In the following sections the actual NASTRAN input is described. The 
section on the Executive Control deck describes the overall system control 
and the available parametric data. The section on Case Control describes 
the control of optional input cases and output requests. The Bulk Data 
section describes the detailed formats for each new bulk data card. The 
final section illustrates through several examples the actual contents of 
the NASTRAN decks for the many optional paths through the program. 
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6.2 HYDROELASTIC EXECUTIVE CONTROL DECK 


The hydroelastic Executive Control deck is similar to that for the 
existing normal modes analysis. Rigid Format 3. When running the hydro- 
elastic analyses, the user must Insert one of the special DMAP ALTER packages 
into his Executive Control deck. These alter cards are delivered with the 
modified NASTRAN system. 

Two special DIAG's have been added for the hydroelastic analysis. 

DIAG 23 - Prints a list of degrees of freedom Including fluid point 
definitions. For each point, an indication is made 
identifying the sets to which it belongs. 

DIAG 24 - Prints the contents of selected displacement sets. For 

each set , a list of all degrees of freedom belonging to the 
set is given. 

These two DIAG's produce similar output to DIAG's 21 and 22 except the 

following hydroelastic sets are included or modified: 

= Structure point 

Uy = Fluid point 

= Free surface point 

U = U + U, 

2 X fr 

= a bits (structure only) 


ab 

= Interior fluid points 

U = U , + U, 
a ab fr 


Hydroelastic DMAP Alters 

Two sets of DMAP alters to Rigid Format 3, Figure 1 and 2, are provided 
to perform the three-dimensional hydroelastic analysis. These alters will 
obtain the hydroelastic solution with either direct or modal formulation. 
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USETFjUSETSjA FjDK GG $ 

f VEC 

1 CHKPNT 

PARTN 

USETF/PV1/»G>«/»X»/»Y» $ 
PV1 s 

KGG*PV1^/KXX,> jKYY $ 

- 

PARTr 
I PARTN 
• EQUIV 

MGGjPVI ji/MXXj, > S 

RG^PVl>yPX/>j /1 s 
RX^RG $ 


CHXPNT 

r PARTN 

1 CONO 

R6 S 

AFj PV1 >/*^AXY jAYY t 
L8L69/NOGRAV $ 


PARTN 

OKGGjPVI^/DKXXjj j DKYY $ 


» CONO 

LBL6<5jNOFREE $ 


I VEC 

USETF/PV2/»Y»/J»FP»/»COHP» S 


CHKPNT 

PV2 % 


* PARTN 

AYYi^PVZ/AFRYjj^/O * 


1 PARTN 

0KYY^PV2:r/DKFRFR^ S 


* LABEL 

LBL69 $ 


CHKPNT 

AFRY^AXYjKYYjDKXXjDKFRF R $ 


t SWITCH 

USET^USETS// $ 


J CHKPNT 

USET»LISETS * 


COMPCFF 

NOSTRUC^OLDSTR S 


1 COMPON 

2iDIFSTIF $ 


I PA RAMP 

//«rOMPL EXk// OIFSCALE=1 •0/0.0/DIF SCAL /// 

$ 

ADD 

KXX>KOGG/KGG//OIF SCAL S 


^ COMPOFF 

I^DIFSTIF S 


I EOUIV 

KXX^KGG S 


* EOUIV 

MXX.MGG S. 

- 

CHKPNT 

KGG^HGG 9 


1 ALTER 

85 


1 LABEL 

NOSTRUC * 


PURGF 

DKAA/NOGRAV S 


1 CONO 

LBL73jN0GRAV $ 


1 EOUIV 

OKXXjDKNN/MPCFI s 


CHKPNT 

OKNN $ 


1 CONO 

L6L71jMrCF2 $ 
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DKNNjDKFF/SINGLE 5 


» CHKPNT 

DKFF $ 
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1 COnO 
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USET>GO>DKFF/DKAA S 


, LABEL 
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1 CHKPNT 

DKA A , . 

— . - - _ 


FIGURE 1, DIRECT FORMULATION DMAP ALTERS 
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FIGURE 1. DIRECT FORMULATION DMAP ALTERS (Cont’d) 
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FIGURE 2. MODAL FORMULATION DMAP ALTERS 
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FIGURE 2. MODAL FORMULATION DMAP ALTERS (Cont'd) 
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Severn! optional parameters may be specified by the user for each type 
of solution. These parameters are used to; (1) control the optional compu- 
tation paths, (2) specify numerical factors to be used in the formulation, 
and (3) allow blocks of DMAP code to be turned "off" for restart from a 
previous checkpoint run. 

These parameters are entered in the Bulk Data deck using the PARAM card 
in addition to those already provided for Rigid Format 3. 

Direct Formulation Parameters: 

1. C0MPTYP - optional - default = -1 

Controls the type of compressibility calculations performed. A negative 
integer will cause a finite compressibility as defined by the KC(SMP para- 
meter. A positive integer will cause a consti iint equation to be generated 
to provide pure incompressibility. 

2* KC0MP - optional - default =1.0 

Value defines the overall compressibility of the fluid volume. The 
definition is fluid bulk modulus divided by total volume. 

3. DIFSTIF - optional default = 1 

A negative integer value causes the differential stiffness matrix to 
be Included for ullage pressure effects. This matrix is availaole from 
the checkpoint tape cf a RJgid Format 4 solution run of the structure 

uiuuel . 

4. DIFSCALE - optional - default = 1.0 

The differential stiffness matrix may be multiplied by the real value 
of this parameter. 

5. NEWM0DE - optional - default = 1 

A negative integer will cause all DMAP statements and alters up to the 
eigenvalue extraction to be skipped. This allows the uuer to restart 
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the original solution to obtain difrent eigenvectors without changing 
the DMAP alter deck. 

6. 0LDSTR - optional - default = 1 

A negative value will cause most structure-related processing to be 
skipped. This allows the user to restart a previous solution, either 
hydro or structure only, and change the fluid model without recomputing 
the unchanged structure. 

Modal Formulation Parameters: 

1. KC0MP - optional - default = 1.0 
(same as direct formulation parameter) 

2. DIFSTIF - optional - default = 1 
(same as direct formulation parameter) 

3. DIFSCALE - optional - default = 1.0 
(same as direct formulation parameter) 

4. NL,vM0DE - optional - default = 1 
(same as direct formulation parameter) 

5. 0LDSTR - optional - default = 1 
(same as direct formulation parameter) 

6. LM0DES - optional - default = -1 

This integer value specifies the number of the lowest structure modes 
to be used when formulating the hydroelastic matrices. A negative value 
indicates all available modes are to be used. 
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6.3 HYDROKLASTIC CASE CONTROL DECK 

The Case Control data for normal modes analysis, Rigid Format 3, is not 
modri-ied for direct hydroelastic solutions. For modal formulation, the data 
is similar except that two sets of subcases must be provided. The first set 
must select the EIGR card, METHOD* , to be used in eigenvalue extraction for 
the structure-only model. Several subcases may be used to define output 
requests for different vectors with the M0DES card. A second set of subcases 
is also needed to define the eigenvalue extraction and output request of the 
combined fluid/structure model. If the NEWM0DE or 0LDSTR parameter is used 
with modal formulation, only the second set of subcases, used for the complete 
model, is required. Several sample Case Control decks used for each formula- 
tion are shown below. 

Direct Formulation: 

TITLE * 

SPC = 10 
METH0D =50 
DISP = ALL 

Modal Formulation; 

TITLE = 

SPC = 10 
SUBCASE 1 

LABEL = M0DES 0F EMPTY STRUCTURE 
METH0D =10 
DISP = N0NE 
SUBCASE 2 

LABEL = MODES WITH FLUID INCLUDED 
MEIHOD =20 
DISP = ALL 
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Modal Fomulation with Selective Output Requests: 

TITl.E “ 

SPC = 10 
SUBCASE 1 

LABEL - STRUCTURE M0DES 1 & 2 
METHOD =» 10 
DISP = ALL 
M0DES = 2 
SUBCASE 3 

LABEL - STRUCTURE M0DES 3 & 4 
DISP = N0NE 
SUBCASE 5 

LABEL - FLUID/STRUCTURE M0DES 1-3 
METH0D =20 
DISP = ALL 
M0DES = 3 
SUBCASE 8 

LABEL - FLUID/STRUCTURE M0DE 4 
DISP = N0NE 

In this example the eigenvectors for only the first two structure modes 
and the first three combined modes will be printed. 

Hydroelastic Output Control 

The structure printout and plotting Case Control requests are used to 
control both the fluid and structure outputs. The following data is 
available : 

1. Structure-related data such as displacements, forces, and stresses 
are processed with normal NASTRAN control. 

2. Fluid internal pressures are output by including their grid point 
Identification numbers in the DISP = output request. If the fluid 
point is on a free surface defined by CFREE data, the actual free 
surface displacements will be printed. 
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3. Both structure and fluid elements may be plotted as undeformed shapes. 
The interior fluid point degrees of freedom are actually pressures 
and should not be plotted as defonrod shapes. 

4. The deformed shape of the free surface may be plotted using the 
"SHAPE" or "VECT0R" plot options. It is recommended that PL0TEL 
elements be used to define Che free surface. If the fluid elements 
CFHEXl, CFHEX2, etc. were used in the requested plot set, all of 
their boundaries are plotted and result in a confused plot. An option 
is available in Che mesh generator program, MESHGEN, to automatically 
generate the free surface PL0TEL data. 

5. The use of modes parameter to control output requests is described 
under the Case Control section. 
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6.4 HYDROELASTIC BULK DATA 


The new Bulk Data cards used for three=»dimensional hydroelastic modes 
analysis are described on the following pages. These cards are used to define 
the fluid and fluid/structure interface. The tank walls and supporting struc- 
ture are defined with NASTRAN structure elements. The actual tank walls must 
be defined by two-dimensional membrane, panel, or plate elements. 

In addition to the new data cards, the following MASTRAN data cards are 
used for special hydroelastic purposes: 

1. GRID cards are used to define the fluid points. Fluid grid points 
contain only one degree of freedom and may not be connected to the 
structure elements. 

2. GRAV cards are used to define the magnitude and direction of the 
gravity field. The set identification numbers are referenced by 
the fluid boundary data cards. 

3. SPC and SPCl data cards may be used to define constraints on the 
fluid grid points. These constraints are used to d' fine regions of 
^ero pressure in the fluid, such as a free surface without gravity 
effects or an anti-symmetric boundary condition on a plane of sym- 
metry. Only degree-of-freedom number 1 may be specified for a 
fluid grid point. 
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BULK DATA DECK 


Input Data Card CFFREE Free Fluid Surface 

Description : Defines those fluid elements composing the free fluid surface 

in a hydroelastic analysis. 


Format and Example : 


1 

^ 2 

3 

^ - 

4 


7 

■V 

8 


10 

CFFREE 

EIDF 

GRAVID 

FACE 


EIDF 

GRAVID 

FACE 

X 

X 

CFFREE 

100 

100 

3 


101 

100 

4 




Field 


Contents 

EIDF 


Fluid element identification number (Integer > 0) 
(see Remark 1) 

GRAVID 


Idfcncif ication number of a GRAV gravity vector set 
(Integer > 0) 

FACE 


Identification number of the face of the fluid element, EIDF 
forming the free surface (0 < Integer 1 6) (see Remark 2) 

Remarks : ' 

1. 

Allowable fluid element types are CFHEXl, CFHEX2, CFTETRA, 
CFWEDGE. 


2. 

The numbering conventions for solid faces are defined in the 
fluid element connection Bulk Data description. 
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BULK DATA DECK 


Input Data Card CFHEXi Fluid Hexahedral Element Connection 


Description : Defines two types of fluid hexahedral elements (three- 

dimensional solid with 8 vertices and 6 quadrilateral faces) to be used in 
hydroelastic analysis. 


Format and Example : 


123456789 10 


CFHEXI 

EID 

MID 

G1 

G2 

G3 

G4 

G5 

G6 

abc 

CFHEX2 

15 

100 

1 

2 

3 

4 

5 

6 

ABC 


+bc 

G7 

G8 








+BC 

7 

8 









Field 


Contents 


CFHEXI 

EID 

MID 

G1 t • • • jG8 


Remarks : 1 . 


CFHEXI or CFHEX2 (see Remark 4) 

Element identification number (Integer > 0) 

Material identification number (Integer > 0) 

Grid point identification numbers of connection points 
(Integer >0, G1 G2 ... ^ G8) 



Element identification numbers must be unique with respect 
to all other element identification numbers. 


2. The numbering and order of the grid points and faces, required 
for specifying free fluid surfaces, are defined in the figure 
above . 


3. The quadrilateral faces must be nearly planar. 

4. CFHEXI is developed by 5 tetrahedra, CFHEX2 by 10 overlapping 
tetrahedra. 

5. Material ID must reference a MAT4 Bulk Data card. 
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BULK DATA DECK 


Input Data Card CFLSTR Fluid /Structure Interface 

Description ; Defines fluid/structure Interfaces for the hydroelastic 
analysis. 


Format and Examples ; 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

CFLSTR 

EIDF 

GRAVID 

EIDSl 

EIDS2 

EIDS3 

EIDS4 

EIDS5 

EIDS6 

abc 

CFLSTR 

100 

10 

1 

2 

11 

12 

21 

?2 

ABC 


+bc 

BIDS 7 

EIDS8 

— 

- etc. - 

— 




dcf 

+BC 

31 

32 







1 


Alternate Form 


CFLSTR 

EIDF 

GRAVID 

EIDl 

"THRU" 

EID2 





CFLSTR 

200 

100 

101 

THRU 

106 






Field 


Contents 

EIDF 


Fluid elements Identification number (Integer > 0) 
(see Remark 3) 

GRAVID 


Identification number of a GRAV gravity vector set 
(Integer > 0) 

EIDSl,, .. 


Structural element identification number (Integer > 0) 

Remarks ; 

1. 

As many continuation cards as desired may appear when "THRU" 
is not used. 


2. 

All element ID's between EIDl and EID2 must exist when using 
"THRU" option. 


3. 

Allowable fluid element t3T>es are CFHEXl, CFHEX2, CFTETRA, and 
CFWEDGE. 
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BULK DATA DECK 


Input Data Card CFTETRA Fluid Tetrahedral Element Connection 

Description ; Defines a fluid tetrahedral element (three-dimensional solid 
with 4 vertices and 4 triangular faces) • to be used in hydroelastic analysis. 

Format and Example ; 


1 23456789 10 


CFTETRA 

EID 

MID 

G1 

G2 

G3 

G4 



>< 

CFTETRA 

25 

100 

1 

2 

3 

4 





Field 

EID 

MID 

G1,G2,G3,G4 


Contents 

Element identification number (Integer > 0) 

Material identification number (Integer > 0) 

Grid point identification numbers of connection points 
(Integer > 0, G1 G2 4 G"i G4) 





Remarks ; 1. Element identification numbers must be unique with respect to 

all other element identification numbers. 

2. The numbering of the grid points and faces, required for 
specifying free fluid surfaces are defined in the figure above. 

3. Material ID must reference a MAT4 Bulk Data card. 
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BULK DATA DECK 


Input Data Card CFWEDGE Fluid Wedge Element Connection 


Description ; Defines a fluid wedge element (three-dimensional solid, with 
three quadrilateral faces and two opposing triangular faces) to be used In 
hydroelastlc analysis. 


Format and Example ; 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

CFWEDGE 

EID 

MID 

G1 

G2 

G3 

G4 

G5 

G6 


CFWEDGE 

25 

100 

1 

2 

3 

4 

5 

6 



Field 


Contents 


EID 

MID 

G1 , • . • , G6 


Element identification number (Integer > 0) 

Material Identification number (Integer > 0) 

Grid point identification numbers of connection points 
(Integers > 0; G1 G2 G6) 



Remarks ; 1. Element identification numbers must be unique with respect 

to all other element identification numbers. 


2. The numbering of the grid points and faces, required for 
specifying free fluid surfaces, are defined in the figure above. 

3. The quadrilateral faces must be nearly planar. 

4. Material ID must reference a MAT4 Bulk Data card. 
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BULK DATA DECK 


Input Data Card MATE Fluid Material Property Definition 

Description : Defines the fluid density for a hydroelastic analysis. 

Format and Example : 


123456789 10 


MATF 

MID 

Q. 

p>< 


1 ^ 





MATF 

103 

.6 




iB 

■Hi 




Field 


Contents 


MID Material identification number (Integer > 0) 

p Mass density (Real > 0.0) 


Remarks : 1. The material identification number may be the same as a MATl , 

MAT2, or MAT3 card, but must be unique with respect to other 
MATF cards. 
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Q Fluid Points 
□ Structure Element 


FIGURE 3. SAMPLE HYDROELASTIC MODEL 



6.5 SAMPLE HYDROELASTIC PROBLEMS 


In this section a set of sample hydroelastic problems are shown which 
demonstrate the various options available to the user. In all the problems, 
the same fluid/structure model Is used. This model, representing a half-model 
of a simple box filled with fluid. Is shown In Figure 3. The structure walls 
are modeled with CQUAD2 elements and the fluid consists of a single CFHEX2 
element . 

Problem I - Direct Formulation Solution 

Figure 4 shows a NASTRAN deck for the solution of the sample hydroelastic 
Problem I.l using the direct formulation method. In this problem, as In all 
following problems, the fluid grids and elements are number 1-8 and the 
structure grids and elements are numbered 101-108. The CFLSTR and CFFREE 
cards are used to define the structure/fluid boundary and the free surface 
respectively. This problem has been checkpointed to allow restarting for 
additional eigenvectors. 

Figure 5 shows the restart of Problem I.l to obtain new modes. The 
NEWM0DE parameter has been added to skip the recnmputatlon of the unchanged 
model. In particular, the NEkM0DE parameter turns off the hydroelastic alters 
which would normally be re-executed during a restart. When the NEWM0DE para- 
meter is used, the only changes that can be made are those which control the 
eigenvalue extraction and the output requests. 

Note that when performing this type of restart, the restart dictionary 
should be stripped back to the start of eigenvalue extraction. For the direct 
formulation DMAP, this occurs at statement 96. This is to assure that the 
proper data blocks from the old problem tape are used In output recovery. 
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I 

I 

I 

I 

I 


r 


I 


I 


I 

I 

1 


TO HYORO^UAl 
TIME 1 

aiAC f 

CHKPNT YES 
SOL S>0 

ARP_ PISP 

$ 

S (DIRECT FORMULATION ALTERS) 

cInd ~ " 

TITLE = HYDRO ELASTIC DIRECT FOR 
„SU5 T_I J L£_^P R OB L£_M_Ij?1_^m_SO. 


MULATION TEST 

Cy Tip N_U IJ H__C H EC KP p I hT 


SPC = 


IQ 


METHOD 
OISP = 
fPCF = 
BEGIN BULK 


= 50 

all„ 

ALL 


GRID 

1 


J3...P 

: .0 

0.0 


GRID 

2 


6*0 

0.0 

0.0 


GRID 

3 


0«Q 

12.0 

0.0 

GRID 

4 


6*0 . 

12.0 

0.0 


GRID 

5 

• 

0*0 

0.0 

12.ti 

- 

GRID 

6 


6<0 

0.0 

12.0 


GRID 

7 



12.0 

12.0 


GRID 

8 


6*0 

12.0 

12.6 


GRID 

101 


0.0 

0.0 

0.0 


GRID 

_102 


6 « 0 

0.0 

0.0 


GRID 

103 


0.0 

12. 0 

oVo 


GRID 

104 


6.0 

12.0 

0.1 


GRID 

-J05 _ 


0.0 

0.0 

1 2.0 


GRID 

106 


6.0 

0.0 

12.0 


GRID 

107 


0.0 

12.0 

12.0 


GR iO 

103 


6.0 

12.0 

12.0 


COUA 02 

101 

100 

101 

10? 

106 

105 

C0UA02 

102 

TOO 

10 2 

104 

108 

106 

C0UAD2 103 

100._ 

104 

103 

107 

108 

C0UA02 

104 

100 

loi 

103 

104 

102 

CFHcXZ 

1 

200 

1 

2 

4 

3 5 

«C1 

a 

7 





CFFREE 

1 

100 

6 




CFLSTR 

1 

100 

101 

THRU 

104 


P0UA02 

100 

100 

_.g,6 




MAT1 

100 

1C.6<6 


.7 

~T92^3" 


MA TF 

200 

9.355-4 





f SPC1 

10 

1256 

101 

103 

105 

__io7 : 

OM IT 1 

4 

101 

10 3 

105 

107 


Q.N I T 1 

4 56 

102, 

104 

106 

108 


GRAV 

100 


386.0 

0.0 

0.0 

-1.0 

EIGR 

50 

GIV 

0.0 ■ 

' 20.0 * 

6 '■ 

6 ■ 0 

«E1 

MAX 









ENODATA 


FIGURE 4. PROBLEM I.l - DIRECT FORMULATION SOLUTION 
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* hS 



10 HYDR-0#UAJL 

e* , . X '' ' 

TIME 1 
0IA6 8 


$ RESTART dTcTIONARY FROM PROBLEM I.l 

$ NOTE •«« OICTIONARY SHOULD BE STRIPPED BACK TO STATEMENT NO* 85 



SOL 3^0 


- 


& 

\ 

APP OISP 




«: 

S 





$ CDIRECT FORMULATION ALTERS) 


1 

* 




il.' 

CENO . 


. 



TITLE = HYDRO ELASTIC DIRECT FORMULATION TEST 


\ 

SUBTITLE = PROBLEM I 

• 2 - RESTART OF SOLUTION FOR 

ADDITIONAL MOOES 

•i 

ECHO both 



* ~ 


SPC = 10 


* 

’ - 


METHOD. = 50 




1 

DISP = ALL 




• 

SPCF ^ ALL 





BEGIN BL'LW 



" 


t 


. 


■* 

$ NEV EI6R CARO 

FOR 

DIFFERENT MODES 

*T* 


$ 




T '' 

/ 9 

10 


— « 

•? 

A 

EIGR 50 

G!V 

100*0 2500*0 

0 SE 


<E1 MAX 





S 




t 

> 

$ PARAMETER TO 

SKIP 

UNNEODEO OMAP 

•rr • 

* 

S 





PARAH NEUMODE 

-1 




ENDDATA 




4 > 

* 



• 

• 



) 

4 





I 

I FIGURE 5. PROBLEM 1.2 - RESTART OF SOLUTION FOR NEW MODES 

\ 
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Proble m 1 1 - Compu tatio n of Separate FLU ID and STRU CTURE Models 

In this problem, the structure model will be computed separately from 
the fluid. This case can occur when it is desired to obtain the behavior 
of an empty tank and then to add fluid at a later time. 

Figure 6 shows the NASTRAN normal modes solution of the structure only. 
Problem II. 1. This run has been checkpointed to allow rest.arting for a 
hydroelastic analysis. 

In Figure 7, the run is shown which will add a fluid model to this 
structure and obtain the solution of the combined model, using the direct 
formulation method. Since the bulk data for the structure model is on the 
old problem tape, only tlie new cards required to define the fluid are provided. 

The 0LDSTR parameter is used to skip rc-c .mputation of the structure. If 
the 0LDSTR parameter is used, no ch.ange to the structure or its constraints 
may be made. 
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ID H YDR0>UAI 
TIHE 1 
CHKPNT = YES 
OlAG 8 

_i^PP_ OISP 

SOL 3i0 
CENO 



II.IL.E.=_HYDROELASTIC _RESTART_ TEST _ . 

SUBTITLE = PROBLEM II. 1 - NORMAL MODES OF STRUCTURE ONLY 
LABEL = NASTRAN SOL 3^0 SOLUTION 
£P C = 10 


METHOD = 50 
OISP = ALL 
-SPCF =_ALL _ 



begin bulk 






J 

GRID 

101 


0.0 

0.0 

0.0 


1 

GRID 

102 


6.0 

0.0 

0.0 



GRID 

103 


0.0 

12.0 

0.0 


X- 

GRID 

1Q4 


6.0 

12.0 

0.0 


GRID 

105 


0.0 

0.0 

12.0 



GR 10 

106 


6.0 

0.0 

12.0 



GRID 

107 


0.0 

12.0 

12.0 


-1- 

tRI_P. 

108 


6.0 

12.0 

12.0 


1 

CD 11 A D2 

101 

100 

101 

102 

106 

105 


COUA D2 

102 

100 

102 

104 

108 

106 

T 

COUA 02 

103 

100 

104 

103 

107 

108 

ar 

COUA D2 

104 

100 

101 

103 

104 

102 


POUA D2 

100 

100 

• 06 




r 

MATl 

100 

10.646 


.3 

.92-3 



SPC1 

10 

1256 

101 

103 

105 

107 


OMIT 1 

4 

101 

103 

105 

107 



OMIT 1 

456 

102 

104 

106 

108 


• 

EIGR 

50 

GIV 

0.0 

20.0 

6 

6 0 


AE1 MAX 

ENDDATa 


«E 


4 


1 


. FIGURE 6. PROBLEfl II. 1 - NOR>»AL MODES OF STRUCTURE ONLY 
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ID HYDROiUAI 
TIME 1 
OIAG B 
$ 

$ RESTART. DICTIONARY FORM PROBLEM II. 1 _ 

$ 

APP OISP 

..SOL-JbO . 

s 

$ CDIRECT FORMULATION ALTERS) 

CENO 

TITLE = HYDROELASTIC RESTART TEST 

-SU8Tj:TLE-.=-.PR0aLEM_XI..2_-_j?£STJlR.T OE_S TURCTURE _ONLY RUN_._ 
LABEL = NASTRAN HYDROELASTIC DIRECT FORMULATION DMAP 
ECHO = BOTH 

SPC_rr-JjO 

METHOD = 60 
DISP = ALL 
SPCF = A LL 


1 BEGIN 

BULK 




4 * 

S KX it 

NOTE - 

S J R »JC T U.R E B.U L K_D A lA 

-.JS-..0N_ 

. R ES T ART_T A P E_„ 

, $ 



" GRID 

1 

0.0 

O.Q 

0.0 

* GJ?I-D. 


6.0 

Q.O 

0.0 

GR ID 

3 

0*0 

12.0 

0.0 

GRID 

4 

6.0 

12.0 

0.0 

GRID 

5 

n.o 

0.0 

12.0 

GRID 

6 

6.0 

O.Q 

12.0 

GRID 

7 

O.Q 

12.0 

12.0 

grid 

8 

6*0 

12.0 

12.0 


CFHEX2 

1 

200 

1 

2 

4 

3 

5 

«C1 

8 

7 






CFf REE 

1 

mo 

6 





CFLSTR 

1 

100 

101 

THRU 

104 



MATF 

200 

9.355-4 






GRAV 100 


386.0 

0.0 

0.0 

-1 .0 


EIGR 

60 

GIV 

0.0 

20.0 

6 

6 

0 

«E2 

MAX 








S parameter TO USE INCOMPRSSIBLE CALCULATIONS 

$ 

PARAK C0MPIYP_1 


S 

$ parameter TO SKIP R F.COMP UTA T 10 N Oh UNCHANGED STRUCTURE 
PARAM OLDSTR -T 

-.END DATA . - 


FIGURE 7. 


PROBLEM II. 2 - RESTART WITH FLUID MODEL ADDED 
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Problem III - Addition of Internal Ullage Pressures 


The NASTRAN differential stiffness analysis can be used to determine the 
effect of an intetnal ullage pressure on a tank filled with fluid. To per- 
form this analysis, two runs are required. The first run is a NAST.^N 
Differential Staffness run on the tank structure only. A sample run. 

Problem III.l, is shown in Figure 8. In this run, a DILAP EXIT instruction 
was altered in to stop computation once the differental stiffness matrix 
was computed. Even though the solution was not obtained, the second subcase 
ill Case Control is still required. The Sol 4,6 subset was used *t> turn off 
all CHKPNT instructions in the Rigid Format. A C'.ni’NT instruction was then 
altered in before the EXIT to save the Differential Stiffness matrix on the 
New Problem tape. This procedure reduces cost since only the one desired 
data block is copied to the new problem tape. Note also that an internal 
pressure of 1.0 was used. This value can be scaled to any d' xred value 
in the second run. 

The hydroelastic analysis solution is shown in Figure 9. In this 
particular run, the direct formulation alters were used, but the modal 
formulation could also be executed. Since the structure bulk data is on 
the restart tape, only the cards needed to define the fluid are used. The 
DIFSTIF parameter is included to trigger addition of the differential stiff- 
ness to the structure. This matrix is obtained from the old problem tape and 
will be scaled by the DIFSCALE parameter prior to addition. 
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I 

I 

I 

I 

1 

1 


I 


ID HYDHOjUAI 
TIME 1 

CHKP^17 = YES 

DIAG a 

APP DISP 

SOL A >6 
ALTER 118 

-CHKPNT KDGG^-S - . - . 

EXIT J 
ENOALTER 

CEHD .. . 

TITLE = HYDROELASTIC UL 
SU3 TITLE = PROBLEM III. 

SPC--= 10 .. 

LOAD = 10 
DISP ALL 

.SUBCASE! 

LAPEL = 

SUBCASE 2 

LABEL-. 

BEGIN BULK 


LAGE PRES 
1 - DIFFE 


SURE TEST 
RENTIAL STIFFNESS 


RUN 


STATIC SOLUTION 


DIFFEPENTlAL.,STiFFNESS SOLUTION 


GR ID 

101 


0.0 

0.0 

O.Q 


GR.ID, 

. 10Z. 


6.0 

.. 0.0 „ 

0.0 


GR ID 

103 


0.0 

12.0 

0.0 


GR ID 

104 


6.0 

12.0 

0.0 


.GR.JQ 

.105 .. 


0.0 

Q.O 

12.0 


GRID 

106 


6.0 

0.0 

12.0 


GR ID 

107 


0.0 

12.0 

12.0 


-.GRID... 

108 


6.0 

12.0 

12.0 


COUA D2 

1 01 

TOO 

101 

102 

106 

105 

COUA D2 

1 02 

100 

102 

104 

108 

106 

COUA D2 _ 

103 

100 

.104 

103 

107 

108 

COUA D2 

104 

100 

101 

103 

104 

102 

P0UAO2 

100 

100 

.06 




.... . MAXI. 

..100 

.._1D.656„ 


.3 

-.92-3 


SPC1 

10 

1256 

101 

1G3 

105 

1D7 

PLOA C2 

10 

1.0 

101 

THRU 

104 



.EHDDATA_ 



I 

1 

1 


FIGURE 8. PROBLEM IIT.i - DIFFERENTIAL STIFFNESS, RUN //I 
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ID HYCROjUAI 

time 1 

DIAG 8 

s 

$ R£START-DLC_T.IONARY T.ROM, PROBLEM. III. l._ ... 

S 

APP DISP 

soL._3:i.a — - 

j 

S (DIRECT FORMULATION ALTERS) 

CEND 

TITLE = HYDROELASTIC ULLAGE PRESSURE TEST 

-SU8Tl_T.LE_=_PRDBLE«_lXl-.2--__hORHA.L_M ODES „RE.S TART WITH HYDROELASTIC DMAP_ 

ECHO = BOTH 
SPC = TO 

METHOD-^_5.0 

DISP = ALL 
SPCF = ALL 

BEiULIL8.UJLJ( 

s 

$ *** note - STRUCTURE BULK DATA IS ON RESTART TAPE 


GRID 
GR ID 
GRID 

1 

2 

3 

0.0 

6.0 

0..Q 

0.0 
0.0 
1 ?.o 

0.3 
0.0 
0*0 

GRID 

4 

6.0 

12*0 

0.0 

GR ID 

5 

0.0 

0.0 

12.0 

-GJUO 

6 

6.0 

0.0 

12.0 

GR ID 

7 

0.0 

12.0 

12.0 

GR ID 

8 

6*0 

12.0 

12.0 


CFHEX? 1 2C0 1 2 4 3 5 6 S 

«C1 3 7 

CFFREE 1 100 6 


C.ELSJ.E_ 

1 

IPO 

101 

THRU _ 

_...1 PA- 

HATF 

200 

9.35 5- 

•4 



0MIT1 

4 

101 

103 

105 

ID? 

r.M T T 1 

._,456 

1 CL2 

. 104 

106 

_-..1.08 - 

GRAV 

100 


386.0 

c.o 

0.0 

EIGR 

50 

GIV 

0*0 

20.0 

6 


-A El i^.AX 

$ 

S PARAMETERS TO TRIGGER ADDITION OF ULLAGE PRESSURE 

PARAM DIFSTIF -1 
PARAM DIFSCALE 14*7 

..END DATA 


FIGURE 9. PROBLEM III. 2 - DIRECT HYDROELASTIC SOLUTION WITH 
DIFFERENTIAl, STIFFNESS ADDED 


Problem IV - Modal Formulation Solution 


Figure 10 shows the NASTRAN deck required for a modal formulation 
solution. Problem IV. 1, of the sample hydroelastic model. Two subcases have 
been provided to control the eigenvalue extraction of the structure only and 
the fluid/structure combination. The bulk data provided is identical to that 
for the direct formulation solution in Problem I.l except for the extra EIGR 
card. This run has been checkpointed for future restarts. 

Figure 11 shows the restart of Problem IV. 1 to obtain new modes. The 
NEWM0DE parameter has been added to skip the recomputation of the unchanged 
model. For this restart, the restart dictionary should be stripped back to 
the start of the eigenvalue extraction. For the modal formulation DMA? this 
occurs after LABEL NEl'/2 in the DMAP alter sequence at statement 119. Note 
that when performing this restart the only changes that can be made are 
those which control the eigenvalue extraction and output requests. 

The input deck for "’roblem IV. 3 is shown in Figure 12. This run per- 
froms a restart of Problem IV. 1 and changes the fluid model. In this case 
the single CFHEX2 element has been replaced by two CFWEDGE elements. The new 
fluid boundary has been given on CFLSTR cards. This type of situation can 
occur when it is desired to know the behavior of a tank with different fluid 
levels. The advantage to the restart technique is the solution to the 
unchanged structure model need not be repeated. The 0LDSTR parameter has 
been included to skip these calculations and cause the structure modes to 
be recovered from the old problem tape. Note also that only one subcase is 
given in Case Control. This is because the extraction of eigenvalues for the 
structure only is to be skipped. Because the fluid Bulk Data cards are not in 
the Rigid Format restart tables, it is also necessary to change at least one 
grid card. In Problem IV. 3, fluid grid 1 ^/as replaced with duplicate grid as 
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shown in Figure 12. This change will then cause the proper modules to be 
re-executed and the new fluid model to be processed. 





0*0 

‘ 6 *' 6 ‘ 

0.0 

0 . 0 ^ 

12.0 





104 


.92-3 I 

IDS 107 

107 

108 

0.0 - 1.0 

10 10 0 «E 

6.60 «E 


SOLUTION 



I 

I 



10 HYORO>UA1 
TIME 1 
01 AG 
S 

.JS.RESTART DICTIONARY FROM PROBLEM IV.1 

$ »n*w NOTE »MH« DICTIONARY SHOULD EE STRIPPED BACK TO LABEL NEW2 
S AT DHAP STATEMENT NO. 119 

SOL 3^0 
APP DISP 

~S < M 0^ A L F O^MU L aTiO N ^ L T EP S ) 

$ 


X^AIDL 

TITLE = HYDRO ELASTIC MODAL FORMULATION TEST 

SUBTITLE = problem IV.2 - RESTART OF SOLUTION FOR NEW MODES 

SPC = 10 

OiSP = ALL 
SUBCASE 2 

L A B E L M OD E S WITH FL UID IN CLUDED 

METHOD = 60 
SPCF = ALL 

BEGIN BULK 

S 

S NEW EIGR CARD FOR DIFFERENT MODE 

_S 

/ 11 12 

EIGR 60 GIV 100.0 2500.0 ' 0 

<E2 MAX 

5 PARAMETER TO TURN OFF UNHEEDED DMAP 

S • " 

"PARAM NEWMCToE -1 
ENDDATA 





4 FIGURE 11. PROBLEM IV.2 - RESTART OF SOLUTION FOR NEW MODES 

j 

I 
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ID h'tl.RCiUAI 
TIME 1 

OIAG »r _ _ _ 

s' 

S restart DICTIONARy FROM PROBLEM IV. 1 

s __ _ 

SOL >>d 

APP nisp 

_s _ _ ^ 

s"iMCDAL FORflULATJON^ALTERS) 

S 

CEND 

TiTL'r= HYDRO Elastic modal formulation test 

subtitle = PROBLEM IV. 3 - RESTART OF SOLUTION WITH NEW FLUID MODEL 
ECHO = BOTH 

sycT^ 10 “ ----- . - 

DISP = ALL 

SUBCASE 2 _ 

iABt:L'= modes' WI TH' FLU 10“ INC LUOcD ' 

METHOD = 60 
SFCF = ALL 

b'egin'bulk 

s 

$ NEW FLUID MODEL 

s' 

/ 1 4 


CFWEDGE 1 200 1 23 5 67 

' CFW£' 0 £E '2 20 C' ‘2 ' 4 3 " 6 8 7 

CFFRFE 1 100 5 2 IDO 5 

CFLSTR 1 101 104 

CFL'STR 2 T 02 10'3 104 

s 


S »»* NOTE »»* AT LEAST ONE GRID MUST BE ALTERED IN TO FORCE 
$ ' RErxfCUTTON CF proper MODULES 

S 

/ 14 

GRID 1 ■ .0 .0 . 0 ' 

s 

$ parameter TO SKIP RECOMP UTATiO N OF UNCHANGED STRUCTURE 

s 

PARAM CLDSTR -1 
ENOOATA 


1 - 

I 

I 

FIGURE 12. PROBLEM IV. 3 - FLUID MODEL CHANGE 
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7.1 AUTO>UTED INPUT GENERATOR - MESHGEN 


MESHGEN Is a special purpose NASTRAN mesh/input generator for axisymmetric 
tank structures. It allows both the finite element idealization of the tank 
shell and its three-dimensional fluid contents for use with the NASTRAN hydro- 
elastic analysis capability. The methodology employed requires minimal 
assumptions to maintain maximum generality for such fluid/structur- geometries. 

MESHGEN is implemented as a stand-alone utility that can be used to 
generate basic models, NASTRAN Bulk Data, and simple structural plots. An 
English language-based control structur*’, MESHLAN, has been chosen to simplify 
use understanding of the input generator processing. The syntax and verbal 
constructs of the language are oriented toward the structural analyst using 
terminology that is familiar. With little practice, many models can be 
generated without reference to a user's manual, thereby further reducing 
model preparation time and potential misunderstandings encountered with com- 
plex documentation. 

Capabilities of MESHGEN 

MESHGEN is capable of performing automatically a subset of mesh and. input 
generation functions that will reduce input preparation time by an order of 
magnitude for large models. The specific tasks performed by MESHGEN are: 

1. Fully automated finite element mesh generation for shells and 
solids of revolution, as well as a special three-dimensional trun- 
cated solid. This capability includes the facility to vary the 
mesh spacing. 

2. Fully automated element definition utilizing any of the NASTRAM 
quadrilateral and triangular plate elements and fluid three- 
dimensional elements. 

3. Complete user control over grid point and element numbering. 
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4. Automated generation of stringers with constant properties reinforcing 
a shell at regularly spaced circumferential or longitudinal stations. 

5. Fartially automated definition of permanent single-point constraints. 

6. Automated generation of structure/fluid interaction data and fluid 
free surface data. 

7. Simple structural plotting capability of both three-dimensional 
models and a two-dimensional pseudo-development. 

These functions will eliminate all or most of the input preparation for 
the following classes of NASTRAN Bulk Data: 

1. GRID cards 

2. CONNECTION cards 

3. PLf)TEL cards for fluid free surfaces 

4. CFFREE and CFLSTR cards 

A feature is provided to eliminate the Inconvenience of generating 
erroneous punched output by affording the jser an opportunity to study the 
automatically generated model to determine its suitability before requesting 
punched output. The NASTRAN Bulk Data images are saved on a ditsk or drum 
file that may be accessed at a later time to punch a deck. 

7.2 FUNCTION.^. DESCRIPTION OF MESHLAN 

This section defines all of the MESHLAN statements that are used to 
generate finite element idealizations and NASTRAN Bulk Data for those geo- 
metries defined in Section III. The basic nomenclature used is described 
below. 

Any word(s) capitalized and underlined are keywords and should appear 
exactly as shown in the description. However, only the first four letters 
of a keyword are required. Examples are: 
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M0DEL 


ELEMENTS 

GRAVITY 

Any statement appearing In braces will be supplied by the user, as in: 

{fluid level value} 

{integer ID} 

{value} 

The use of square brackets denotes a set of options that may be selected by 
the user, i.e. : 

TANK 

FLUID 

GFLUID 

TFULL 

R 

THETA 

Z 

Certain orderings of MESHLAN st:- cements must be followed to avoid ambiguity. 
All statements are defined by a hierarchical level. Commands at the same 
level may be ordered in any manner, but subcommands must appear beneath their 
higher level command. A brief outline of the hierarchical structure is shown 
below (the numbers in ( ) refer to the number in the command description 
appearing later) . 
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M0DEL (1) 


GC0METRY (2) 

SHAPE (3) 

BOUNDARY (4) 

MESH (5) 

0 PR0PERrY (6) 

STEP or DIVIDE (7,8) 

ZPR0PERTY (9) 

SHELL or S0LID (1C, 11) 

NUMBER (12) 

INSYS (13) 

0UTSYS (14) 

FIX (15) 

ELEMENTS (16) 

PR0PERT Y (17) 
THICKNESS VARIES (18) 

STRI NGERS (19) 

AL0NG (20) 

GRiWITY (21) 

PL0T (22) 

PLT2 (23) 

(23n) 

RfNtm (24) 

F IND B0UNDARIES (25) 



1. LEVEL 1 COMMAND M0DEL 


SYNTAX: 

M0DEL {model name} [ > aVE [.NEW] ) 

Defines the model Identification, (model name} may be any string of 8 
or less BCD characters starting with a letter. The SAVE option causes 
Bulk Data images to be written to the SAVE disk file for deferred 10 
operations. (See Section V for a detailed description.) The first 
time a model is saved on a given disk file, NEW must be specified. NEW 
may not appear without SAVE . 

EXAMPLES: M0DEL TANKlOl 

M0DEL SRU, SAVE, NEW 
M0DEL SHUTTLE, SAVE 

NOTE: Every case must begin with a M0DEL command. 

2. LEVEL 1 COMMAND GE0METRY 
SYNTAX; 


GE0METRY 


Defines the start of a sequence of geometry specification commands. 
NOTE: This command must be present if a model is being generated. 


3. LEVEL 2 COMMAND SHAPE 


SYNTAX; 


SHAPE 


SHELL 0F REV0LUTI0N 


TANK 

S0LID 0F REV0LUTI0N 


FLUID 

TRUNCATED S0LID 


GFLUID 

C0NTAINED S0LID 

OR 

TFULL 

CAPSHELL 


CAPS 

CAP FLUID 


CAPF 

CAPB0TH 

L 


CAPB 


5 





I 


I 

I 


4 


? 


« -‘ 


i 


I 

I 

I 

I 


\ 

’ .. 


Defines the basic geometric configuration of the model. (See Section III 
for a detailed description of the basic shapes.) 

NOTE: This command may only appear as a subcommand of G E0METRY . 


4. LEVEL 2 COMMAND B0UNDARY 
SYNTAX: 


B 0UNDARY 


FUNCTI0N (function ID} 

TABLE {table ID) [INTERP0LATE] 


[({fluid level value})] 


Defines the boundary of the body of revolution either as an explicit 
function or a table. I'/hen the user specifies an explicit function, he 
may select any mesh divisions for his model. If TABLE is specified, 
and INTERP0LATE is absent, the alues of Z given in the table will 
define the longitudinal mesh lines. If INTERP0LATE is present, then any 
divisions may be requested. The {fluid level value} is used to specify 
the axial position of the free surface of the fluid in the tank. This 
value Is required to generate PL0IEL and CFFREE Bulk Data for any SHAPE 
involving a fluid, and must coincide with an axial station. 

EXAMPLES : 

B0UNDARY = TABLE 2 
B0UNDARY » FirNCTI0N 12 (3.0) 

B0UNDARY = TABLE 8 INTERP0LATE (5.5) 

NOTE: This command may only appear as a subcommand of GE0METRY. 

A GRAVITY command (number 21) must be present if a fluid level 
is given. 

TABLE without INTERPOLATE may not be used for SHAPE =• GFLUID. 
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5. LEVEL 1 COMMAND MESH 


SYNTAX; 

MESH ({[Gl] [,G2]}) 

Defines the initial grid point ID's for the raodel and signifies the 
start of the mesh definition commands. Gl is the initial ID for the 
shell portion of the model, G2 for the fluid. Gl, G2 or both will be 
present as required by the SHAPE . 

EXAMPLES : 


MESH (101) 

MESH (1,101) 

MESH (,1001) 

NOTE; A MESH command must be present if a model is being generated. 

6. LEVEL 2 COMMAND 0PR0PERTY 
SYNTAX: 


0PR0PERTY = {id} 

Defines an overall property for the entire model. 


7. LEVEL 2 COMMAND STEP 
SYNTAX: 


STEP 


R 

THETA 

Z 


FR0M {value 1} T0 {value 2} ^ {increment} 


Defines a zone of the model bounded in the specified direction by 
{value 1} and {value 2}. The integer {increment} defines the number 
of elements across the zone. 
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NOTE: STEP may appear only as a subcommand of MESH . 


8. LEVEL 2 COMMAND DIVIDE 
SYNTAX: 


DIVIDE 


R 

THETA 

Z 


[ BY {increment}] 


DIVIDE commands perform the same function as the STEP command, only they 
define divisions over the entire coordinate range specified on a boundary 
data card. (These data cards are defined at the end of this section.) In the 
THETA direction, a DIVIDE command assumes a range of 0® S 6 S 360®. The 
{increment} is not required in the Z direction if a TABLE (no INTERP0LATE) 
boundary is specified. 

NOTE: May appear only as a subcommand of MESH . 


9. LEVEL 3 COMMAND 2PR0PERTY 
SYNTAX: 

ZPR0PERTY = {id} 

Defines a zone element property identification for a particular STEP or 
DIVIDE command. The property ID throughout the zone will be taken as: 

PID = 0PR0PERTY + ZPR0PERTY 

NOTE: If 0PR0PERTY and ZPR0PERTY are both absent, the Property ID’s are 
not punched. NASTRAN will assume a Property ID the same as the 


Element ID. 
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LEVEL 3 COMMAND SHELL 


SYNTAX: 

SHELL 

Used only for TFULL or CAPB^TH SHAPES to allow separate attributes to 
be specified for the shell portion of the model. 

11. LEVEL 3 COMMAND S^LID 
SYNTAX: 

S0LID 

Used only for TFULL or CAPB0TH SHAPES to allow separate attributes to 
be specified for the fluid portion of the model. 

12. LEVEL 4 C0MMAND NUMBER 
SYNTAX: 

NUMBER GRIDS BY {increment}, ELEMENTS BY {increment} 

Defines the numbering Increments desired for grid points and elements 
in the zone. 

NOTE: Extreme care must be taken to assure that zone numberings do not 

overlap causing non-unique ID's. 

13. LEVEL 4 COMMAND INSYS 
SYNTAX: 

INSYS - {id} 

Defines an input coordinate system identification number (location) for 
all grid points in the zone. 
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lA. LEVEL A COMMAND 0UTSYS 


SYNTAX: 

0UTSYS » {id} 

Defines an output coordinate system identification (displacements) for 
all generated grid points in the zone. 

15. LEVEL A COMMAND FIX 
SYNTAX: 

FIX (dof code} ^ [{coordinate list}) 

Defines permanent single-point constraints along grid lines within a 
zone. The {dof code} is the normal NASTRAN code of digits 1-6 with no 
embedded blanks. The {coordinate list} is any list of 25 or less values 
of coordinates along which the dof will be fixed. 

EXAMPLE: 

FIX 123A56 AT (3.0) 

FIX 2A6 AT (0.0, 3.0, 6.0) 

NOTE: Permanent constraints may not be applied to fluid points. 

16. LEVEL A COMMAND ELEMENTS 
SYNTAX: 

ELEMENTS = {element type}, {id} 

Defines the element types desired for the zone. The types must be con- 
sistent with the SHAPE of the model, {id} is the initial element ID 


number. Legal element types for each shape are summarized below. 



Elements 


Shape 

TANK QDMEM, QDMEMl, QDMEM2, QDFLT, QUADl, QUAD2, SHEAR, TWIST, 
T8MEM, TRBSC, TRPLT, TRIAl, TRIA2 

FLUID FHEXl, FHEX2 

GFLUID FHEXl, FHEX2 

CAPS TRMEM, TRBSC, TRPLT, IRIAl, TRIA2 

CAPF FWEDGE, FHEXl, FHEX2 

NOTE: ELEMENTS commands may only appear In Z zone definitions. 

17. LEVEL 5 COMMAND PR0PERTY 
SYNTAX: 

PROPERTY - {id} [, {alpha}] 

Defines element properties. The final property ID will be: 

PID - 0PR0PERTY + ZPR0PERTY + PR0PERTY 
{alpha} is the material orientation angle for relative plate elements. 

18. LEVEL 5 COMMAND THICKNESS 
SYNTAX: 

THICKNESS VARIES 

Used to direct the program to use different Property ID*s at each 
coordinate station to allow for uniform thickness variations in plate 
elements . 

19. LEVEL 3 COMMAND STRINGER 
SYNTAX: 


STRINGER 


BAR 

R0D 


{property id}, {id l}, {id 2} 
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Defines STRINGER (R0DS or BAR) that stiffen the shell, {property id} 
must he 0 for BARs and a BAR0R card input by user into the NASTRAN Bulk 
Data deck, (id 1} and {ii 2} are the initial element ID's in the THETA 
and Z directions! respectively. Both numbers must be present. 

EXAMPLES : 


STRINGER - BAR 0, 101, 201 
STRINGER = R0D 100, 1001, 2001 

NOTE: Only constant property stringers are allowed. It is possible to 

vary stringer properties by manually punching the remaining fields 
on the CBAR cards. 


20. LEVEL A COMMAND AL0NG 


SYNTAX: 


AL0NG 


THETA 

Z 


AT ({list of 


Z 

THETA 


coordinate values}) 


Defines the mesh lines along which the stringers run. The {list of 
coordinate valurts? must correspond to grid lines. The stringers will 
extend to the boundaries of the mesh and will be spaced as defined by 
the list. 


EXAMPLES : 


AL0NG Z AT (90.0, 30.0) 
AL0NG THETA AT (1.0, 0.0) 


21. LEVEL 1 COMMAND GRAVITY 
SYNTAX: 

GRAVITY - {id} 

Defines a GRAV ID to be used on free surface (CFFREE) and fluid/structure 
(CFLSTR) Bulk Data cards. 
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22. LEVEL 1 COMMAND PL0T 
SYNTAX: 

PL0T [{a}, {B}][. N0NUM3 

Requests three-dimensional plot of the currant model, {a} and { 3 } 
are angles to change the plot orientation. The natural view angle 
is looking toward the structure at 6 0°, a controls rotations about 

the Z axis, and 3 the Inclination angle as shown below: 

Z 



If N0NUM is elected, no grid ID numbers will be printed. 

NOTES: Be sure appropriate plot tape is mounted. The order of 

rotations Is v. v/llowed by 3* If cx and 3 are absent, they are 
taken as 0.0. 

23. LEVEL 1 COMMAND PLT2 
SYNTAX: 

PLT2 [. N0NUM 3 

Generates a two-dimensional pseudo-developed plot of the current model. 
N0NUM is the same as in PL0T . PLT2 may not be used except for shell of 
revolution and shell cap models. 
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— ADDENDUM -- 


23a. LEVEL 1 COMMAND PLTHEAD 
SYNTAX; 

PLTHEAD NAME, ACC'T,J0BNAME, BIN# 

Used to generate SC4020 plot identification frame. The input values 
are: programmer name, account ntimber, job name, and bin number. 

All values m<‘St be 8 or less BCD characters beginning with a letter. 


NOTE: This command is required to generate plots at MSEC on the 

IBM 360/65. 
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LEVEL 1 COMMAND PUNCH 
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24. 

SYNTAX: 

PUNCH [{model name}] 

Causes the Bulk Data images to be punched on cards. Punching the 
current model does not require {model name). It Is required to retrieve 
a model from the SAVE file. 

25. LEVEL 1 COMMAND FIND B0UNDARIES 
SYNTAX: 

FIND BOUNDARIES {name 1}, {name 2} 

This command causes fluid/structure Interface data (CFLSTR) to be 
generated when the structure and fluid models have been generated 
separately, {name 1} and {name 2} are, respectively, the model names 
of the STRUCTURE and FLUID In that order . Each must exist on the SAVE 
file. Neither model can be of the TFULL or CAPB0TH types. 

26. LEVEL 1 COMMAND ENDM 
SYNTAX; 

ENDM 

Ends the model case. 

NOTE; Must be present for every case run, whether for model generation 
or 10 operations. 

Specifying the Boundary Function 

Tne boundary function defining the shell or solid of revolution may be 
specified either as a table of (r,z) values, or as the incomplete conic: 

+ a2Z + a^r + a^r ■ a^ + a^ 
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The input format for this function definition is free format, i.e., the 
individual values on the card may be placed in any position on the card as 
long as they are separated by a comma or blanks and maintain the order 
defined below (no data may appear beyond column 72). 

1. Data Delimiter 

$DATA 

2. Functional Definition 


Input No. 
1 
2 

3-4 

5 

6-11 


Contents 

FUNCTION 

ID 

2 ^, Z 2 closed Interval defining boundary 
a used for truncated s'‘lld 
coefficients of conic 


3. Tabular Definition 


Input No. 


1 

2 

3-4 


5 

6 


7-8 

9-10 

11-12 

etc. 


Contents 

TABLE 

ID 

* 1 ’ *2 
a 

NPT number of points in table 

function definition (as many cards as 
necessary) . The values and Z 2 must correspond 
to the exact limits of the structure. 
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Control Sequences 





The commands required for a given case are determined by the operations 
being performed. These operations are model generation or 10. The allow- 


able processes are shown below, indicating the required primary conmtands. 

1. Model Generation 

M0DEL 

GE0METRY 


IffiSH 


PL0T 

PUNCH 

ENDM 


optional 


2. Punching Bulk Data from the SAVE Pile 
M0DEL 

PUNCH {model name} 

ENDM 


3. Generation of CFLSTR cards when the shell and fluid models have been 

generated separately and saved. 

M0DEL 

GRAVITY 

FIND B0UNDARIES (shell name}, (fluid name} 

PUNCH 

ENDM 


File Definitions 

MESHGEN requires three or four scratch files, depending on the operation 
being performed, and two output files. The F0RTRAN unit numbers have been 
assigned as shown below, but may be changed by updating subroutine MESHBD. 
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File Always Required ; 


//FT15F001 DD UNIT-SYSDA, SPACE- (TRK, (1,1)) . 

It DCB- (REC FM-FB , LRECL-80 , BLKSl ZE-1 2 9 60) 

//FT19F001 same as above 

//FT18F001 DD UNIT-SYSDA, SPACE-(TRK, (1,1)) , 

// DCB- (RECFM-UBS ,LRECL-84 ,BLKSIZE-13024) 


Required to Perform SAVE : 

//FT17F001 DD DSN-anyname, V0L-SER-id, 

// DISP-(,CATLG), SPACE- (TRK, (a, b)), 

// DCB-(RECFM-VBS,LRECL-84,BLKSIZE-13024), UNIT-3330 

a and b are determined by the amount of data to be saved. For each 

track allocated, 150 Hulk Data images may be saved. 

Required to Perform PLOTS ; 

//FT14FOOI DD UNIT-SYSDA, SPACE-(TRK, (1,1)) , 

/ / DCB- (RECFM-VBS , LRECL-84 , BLKSIZE-13024 ) 

and the SC4020 plot output unit definition such as 

//FT16F001 DD IINIT=(TAPE7) ,LABEL-(1,BLP) , 

// DISP- (NEW, KEEP) ,DCB-(RECFM-F,LRECL-1152, 

// BLKSIZE-1152) ,V0L- (PRIVATE, ,SER-SC4020) 

The SC4020 file is system dependent and varies from installation to 


installation. 
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7.3 >ffiSHGEN GEOMETRY TYPES 

MESHGEN can generate models for seven basic geometric configurations. 

They are; 

1. Shell of Revolution 

A shell is defined by rotating the function r=f(z), defined on the closed 
interval [a,b], through an angle of 6 degrees about the z axis. The resultant 
surface may be open, closed in 6 (i.e., 0 ■= 360“), closed in z (i.e., f(a) = 0 
or f(b) = 0) or closed in 0 and z. Sample shapes are; 
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2 . Solid of Revolution 

A solid is generated by the same procedure as the shell. 

3. Truncated Solid of Revolution 

This figure is defined as an ordinary shell of revolution, but has an 
additional parameter, a. This is the angle the top surface of the solid 
makes with the line z = a, i.e., 

z 



This model is used to define fluids with tilted surfaces. 
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4. Contained Solid 


This option allows a tank (shell of revolution) and Its fluid contents 
(solid of revolution) to be generated in one pass. 

5. Shell Cap 

The shell cap Is a two-dimensional model that corresponds to the shell 
of revolution. The only difference is that the elements are generated in 
an Isoparametric manner. In particular « the axial (z) coordinate values 
are not defined in equal increments, but the actual arc length of the gener- 
ating function is divided equally. Depending on the function, this may allow 
better geometry of the finite elements. For example: 

Shell of Revolution Model: 



Shell Cap Model: 



Note that a shell cap can NOT be closed In 0 and MUST be closed in z. 

6. Solid Cap 

The solid cap Is analogous to a three-dimensional shell cap. It also 
provides for better element behavior in certain cases. This model has an 
isoparametric style cross-section and equal Increments in 0. For example: 
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Solid of Revolution: 


z 


Solid Cap: 




The solid cap MUST be closed in z. 

7. Cap Both 

This option allows a tank and fluid model to be generated in one step 
using the Solid Cap and Shell Cap methods. 

These seven capabilities allow for well-behaved finite element models 
of tank/fluid systems. 


7.4 DEMONSTRATION PROBLEMS AND USER’S GUIDE 

The following section is a collection of sample problems demonstrating 
the MESHGEN capabilities. There are both simple and more complex examples 
to help the user become familiar with this powerful tool. The problems are 
numbered to correspond to the seven geom.etric configurations that may be 
modeled. In addition, the last problem is a deck that generated a model 
actually used to solve the SRI tank problem. 
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SAMPLE PROBLEM 1.1 


1.1.1 Description ; 

Model a 90° segment of a right circular cylinder, with 5 inch 
radius, on the Interval z»[3.,0.] . Assume homogenious prop- 
erties and use QDMEM elements. 


1.1.2 MESHLAN Program ; 

M0DEL TESTl.l 
GE0METRY 

SHAPE = TANK 
B0UNDARY = FUNCTION 1 
MESH(l) 

0PR0P = 100 
DIVIDE Z BY 3 

NUMBER GRIDS BY 5, ELEMENTS BY 4 
ELEMENTS = QDMEM,1 
STEP THETA FR0M 0.0 T0 90.0 BY 4 
NUMBER GRIDS BY 1, ELEMENTS BY 1 
PL0T -45.0,20.0 
PLT2 
ENDM 
$DATA 

FUNCTION 1,3.,0.,0.,0.,0.,0.,1.,0.,5. 


NOTES; a. The DIVIDE Z BY 3 command means that the range of Z 
specified by the FUNCTION card will be used, 
b. A three dimensional and two dimensional pseudo- 
developed plot will be made. 

1.1.3 Model ; 

Figures 1.1 and 1.2 are the MESHGEN plots generated by this problem. 
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SAMPLE PROBLEM 1.2 


1.2.1 Description ; 

2 2 

Generate a model of the upper hemisphere of the sphere z >6z+r 0 

on the Interval z-[3.,6.]. Model the structure with QUAD2 elements 
assuming a uniform axial variation In shell thickness. Also, 
assume that output displacements are requested In another coordinate 
system, CID ■ 200. 


1.2.2 MESHLAN Program ; 

M0DEL TEST1.2 
GE0METRY 

SHAPE = SHELL 0F REV0LUTI0N 
B0UNDARY - FUNCTI0N 2 
MESH(l) 

0PR0P = 100 
DIVIDE Z BY 3 
0UTSYS - 200 

NUMBER GRIDS BY 100, ELEMENTS BY 100 
ELEMENTS = QUAD2,1 
THICKNESS VARIES 
DIVIDE THETA BY 8 

NUMBER GRIDS BY 1, ELEMENTS BY 1 
PLOT 0.,20. 

PLT2 

ENDM 

$DATA 

FUNCTION 2,3. ,6. ,0. ,1. ,-6. ,1. ,0. ,0. ,0. 


1.2.3 Mode^ . 

The MESHGEN plots for this model are shown in figures 1.3 and 1.4. 
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SAMPLE PROBLEM 1.3 



1.3.1 Description ! 

Model a 90° segment of a right circular cylinder r»4 with homogeneous 
properties. In addition, fix all degrees of freedom at z*2 and fix 
dof 246 at 6«0° and 90°. Use TRIAI elements. 



1.3.2 MESHLAN Program ; 


I 

I 

r 

< 

I 


j 

■V. 


M^DEL TEST1.3 
GE0METRY 

SHAPE - TANK 

B0UNDARY - TABLE 48 1NTERP0LATE 
MESH(lOl) 

0PR0P - 100 
DIVIDE Z BY 3 

NUMBER GRIDS BY 1, ELEMENTS BY 8 
ELEMENTS - TRIAI, 1 
FIX 123456 AT(2.0) 

STEP THETA FROM 90.0 TO 0.0 BY 4 

NUMBER GRIDS BY 100, ELEMENTS BY 1 
FIX 246 AT (0.0,90.0) 

PL0T -45., 20. 

PLT2 

ENDM 

$DATA 

TABLE 48,0.,2.,0.,2,4.,0.,4.,2. 


1.3.3 Model : 

Figures 1.5 and 1.6 are MESHGEN plots of this model. 
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SA>iPLE PROBLEM 1.4 
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1.4.1 Pescrlption ; 

2 

Model the 360* surface generated by the parabola 25z*r using TRIAl 
elements and assuming homogeneous properties. 


F 

i 



1.4.2 MESHLAN Program ; 

M0DEL TEST1.4 
GE0METRY 

SHAPE = TANK 
B0UNDARY = FUNCXI0N 3 
MESH(l) 

0PR0P = 100 
DIVIDE Z BY 3 

NUMBER GRIDS BY 8, ELEMENTS BY 1 
ELEMENTS = TRIAl, 1 
DIVIDE THETA BY 8 

NUMBER GRIDS BY 1, ELEMENTS BY 3 
PL0T 0,,20. 

PLT2 

ENDM 

$DATA 

FUNCTI0N 3,1. ,0. ,0. ,0. ,1. ,-4. ,0. ,0. ,0. 


r 1.4.3 Model: 

i 

A 

The model is shown in Figures 1.7 and 1.8. 
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SAMPLE PROBLEM 1.5 


1.5.1 Description ; 

ffodel a 120® segment of a right circular cylinder r*5 on z®[5.,3.] 
with mesh refinements in both the radial and axial coordinate 
directions. Properties vary linearly In axial direction and the 
shell is stiffened by longerons at z«5 and 9*0®, 90® and 120®. 

Use QDPLT elements. 


1.5.2 MESHLAN Program 

M0DEL TEST1.5 
GE0METRY 

SHAPE = SHELL 0F REV0LUTI0N 
B0UNDARY - FUNCTI0N 8 
MESH(l) 

STEP Z FR0M 5.0 TO 3.0 BY 2 

NU^fflER GRIDS BY 7, ELEMENTS BY 50 
ELEMENTS - QDPLT, 1 
THICKNESS VARIES 
STEP Z FROM 3.0 TO 1.0 BY 3 

NUMBER GRIDS BY 7, ELEMENTS BY 50 
ELEMENTS = QUAD2.151 

PROPERTY * 1000,15.0 
STEP THETA FROM 0.0 TO 90.0 BY 3 
NUMBER GRIDS BY 1, ELEMENTS BY 1 
STRINGERS»BAR 0,10001,20001 
AL0NG Z AT (0.0,90.0,120.) 

AL0NG THETA AT (5.0) 

PL0T -60., 20. 

ENDM 

$DATA 

FUNCTION 8,5. ,1. ,0. ,0. ,0. ,0. ,1. ,0. ,5. 


1.5.3 Model ; 

A three-dimensional MESHGEN plot of the model i shown in Figure 1.9. 
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SAMPLE PROBLEM 1.6 


1.6.1 Description ; 

The problem solved in 1.5 is re-executed to Illustrate the capability 
of generating overlapping eleinents of different types in different 
axial zones . 


1.6.2 MESHLAN Program ; 

M0DEL TEST1.6 
GE0METRY 

SHAPE = SHELL 0F REV0LUTI0N 
B0UNDARY = FUNCTI0N 8 
MESH(l) 

0PR0P - 100 

STEP Z FR0M 5.0 T0 3.0 BY 2 
ZP10P * 10 

NUMBER GRIDS BY 7, ELEMENTS BY 50 
ELEMENTS » QDPLT,1 
THICKNESS VARIES 
ELEMENTS » SHEAR, 1001 
PR0PERTY = 50 

STEP Z FROM 3.0 T0 1.0 BY 3 
ZPR0P =20 

NUMBER GRIDS BY 7, ELEMENTS BY 50 
ELEMENTS = QUAD2.151 

PR0PERTY = 1000,15.0 
ELEMENTS = TWIST, 1151 
PR0PERTY = 800 
ELEMENTS = QDMEM2,2151 
PR0PERTY = 655,25. 
THICKNESS VARIES 
STEP THETA FR0M 0.0 T0 90.0 BY 3 
NUMBER GRIDS BY 1, ELEMENTS BY 1 
STEP THETA FR0M 90.0 T0 120.0 BY 3 
NUMBER GRIDS BY 1, ELEMENTS BY 1 
STRINGERS=BAR 0,10001,20001 
AL0NG Z AT (0.0,90.0,120.) 

AL0NG THETA AT 0-0) 

ENDM 


1.6.3 Model ; 

The model resulting, with respect to grid point and element configura- 
tion, is the same as 1.5. 
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SAMPLE PROBLEM 1.7 


1.7.1 Description ; 

Model of a quadrant of a hemisphere where the axial stations are 
explicitly defined by the user. 


1.7.2 MESHLAN Program ; 

M0DEL TEST1.7 
GE0METRY 

SHAPE - TANK 
B0UNDARY = TABLE 88 
MESH (101) 

DIVIDE Z 

NUMBER GRIDS BY 1, ELEMENTS BY 1 
ELEMENTS = QDMEM2.1 
STEP THETA FROM 0.0 TO 90.0 BY 4 

NUMBER GRIDS BY 100, ELEMENTS BY 6 
PL0T -45., 20. 

PLT2 

ENDM 

$DATA 

TABLE 88, 6., 12. ,0.,7,6. ,6. ,5.657,8. ,4.472,10. ,3.317, 
11., 2., 11. 657,1., 11. 916,0., 12. 


1.7.3 Model ; 

MESHGEN plots are shown In Figures 1.10 and 1.11. 
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SA ?LE PROBLEM 1.8 



1.8.1 Description ; 

A model is generated for a 90° segment of a complex tank composed 
of three distinct sections. These are executed as three cases which 
will be merged to generate the final Bulk Data deck. 


1.8.2 MESHLAN Program ; 

M0DEL TEST1.8A 
GE0METRY 

SHAPE - TANK 
BOUNDARY - TABLE 6 
MESH(lOl) 

DIVIDE Z 

FIX 123456 AT (15.) 

NUMBER GRIDS BY 100, ELEMENTS BY 10 
ELEMENTS ■= QUAD2.1 
PR0PERTY =« 100 

THICKNESS VARIES 
STEP THETA FROM 0.0 T0 90.0 BY 4 
FIX 246 AT (0.0,90.) 

NUMBER GRIDS BY 10, ELEMENTS BY 1 
STRINGER = BAR 0,1001,2001 
AL0NG A AT (45.0) 

AL0NG THETA AT (10.0) 

PL0T -45.0,20. 

ENDM 

M0DEL TEST1.8B 
GE0METRY 

SHAPE «=■ TANK 
B0UNDARY “ FUNCTI0N 12 
MESH(601) 

0PR0P » ICO 
DIVIDE Z BY 3 

NUMBER GRIDS BY 100, ELEMENTS BY 10 
ELEMENTS - QUAD2,51 
STEP THETA FR0M 0.0 T0 90.0 BY 4 
FIX 246 AT (0.0,90.) 

NUMBER GRIDS BY 10, ELEMENTS BY 1 
STRINGER = BAR 0,1101,2101 
AL0NG Z AT (22.5,67.5) 

AL0NG THETA AT (3.0) 

ENDM 

M0DEL TEST1.8C 
GE0METRY 

SHAPE - TANK 

B0UNDARY - TABLE 73 INTERP0LATE 
MESH(901) 

0PR0P = 100 


I 


40 








f 



STEP Z FR0M 3.0 T0 0.0 BY 3 

NUMBER GRIDS BY 100, ELEMENTS BY 10 
ELEMENTS - QUAD2.81 
THICKNESS VARIES 
STEP THETA FR0M 0.0 T0 90.0 BY 4 
FIX 246 AT (0.0,90.) 

NUMBER GRIDS BY 10, ELEMENTS BY 1 
PL0T -45., 20. 

PLT2 

ENDM 

$DATA 

TABLE 6, 15., 10., 0.0, 6, 6., 10. ,5. 4, 12., 4. 55, 13. 5, 4.0, 

14.2.3.55.14.6.3.0. 15.0 

FUNCTI0N 12,10.0,3.0,0.0,0.,0. ,0.,1.,0.,6. 

TABLE 73, 3. 0,0. 0,0. 0,7, 0.0, 0.0, 1.0,0. 3,2.0,. 7, 3.0, 

1.1.4.0. 1.6.5.0.2.25.6.0.3.0 


1.8.3 Model ; 

Figure 1.12 shows the boundary definition of the tank structure. 
Three-dimensional plots of the A and C portion, and a two-dimensional 
plot of the C portion are shown in Figures 1.13-1.15. 
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SAMPLE PROBLEM 2.1 

2.1.1 Description ; 

Model the three-dimensional solid segment of the right circular 
cylinder defined in 1.1 using FHF.X2 elements. 


2.1.2 MESHLAN Program ; 

M0DEL THST2.1 
GE0METRY 

SHAPE = FLUID 

B0UNDARY = FUNGI I0N 1(3.0) 

GRAV = 100 
MESH(,1) 

0PR0P =200 
DIVIDE Z BY 3 

NUMBER GRIDS BY 5, ELEMENTS BY 4 
ELEMENTS = FHEX2,1 
STEP THETA FROM 0.0 T0 90.0 BY 4 

NUMBER GRIDS BY 1, ELEMENTS BY 1 
DIVIDE R BY 3 

NU>ffiER GRIDS BY 100, ELEMENTS BY 100 
PL0r -4O.,2O.O,N0NUM 
ENDM 
$DA1A 

FUNCTION 1,3.,0.,0.,0.,0.,0.,1.,0.,5. 


NOTES: a. A free surface value (3.0) and GRAV command are given 

in order to generate PL0TEL and CFFREE Bulk Data, 
b. N0LTJM is specified on the PL0T card to omit grid point 
numbering. 

2.1.3 Moiel : 

Figure 2.1 shows the three-dimensional MESHGEN plot for this solid. 
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SAMPLE PROBLEM 2.2 
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2.2.1 Deacription : 

^ 2 2 
Model the solid of revolution generated by the sphere z -6z + r =0 

on the interval [3., 6.] as in 1.2. Use HEXl elements and specify an 

output coordinate system for displacements CID = 200 


2.2.2 MESHLAN Program ; 


Mf'DEL TESTE. 2 
%E0METRY 

SHAPE = S0LID 0F REV0LUTI0N 
B0UNDARY = FUNCTI0N 2 
MESH(,1) 

0PR0P =200 
DIVIDE Z BY 3 
0UTSYS=2OO 

NUMBER GRIDS BY 100,ELTHENTS BY 100 
ELEMENTS = FHEX1,1 
DIVIDE THETA BY 8 

NUMBER GRIDS BY 1, ELEMENTS EY 1 
DIVIDE R BY 3 

NUMBER GRIDS BY lO.EI.EMENTS BY 10 
PL0T 5.,2O.,N0NUM 
ENDM 
$DATA 

FUNCTI0N 2,3.,6.,0.,1.,-6.,1.,0.»0.,0. 


2.2.3 Model : 

The three-dimensional plot for the model is shown in Figure 2.2. 
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30 VIEW ANCLES 5.0 20.0 OF TESTa.a 


FIGURE 2.2 
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SAMPLE PROBLEM 2.3 


2.3.1 Description : 

Generate a solid model for the function defined in 1.3. 


2.3.2 MESHGEN Program : 

M0DEL TEST2.3 
GE0METRY 

SHAPE = FLUID 

B0UNDARY = TABLE 48 INTERP0LATE 
MESH(.lOl) 

0PR0P =200 
DIVIDE Z BY 3 

NUMBER GRIDS BY 1, ELEMENTS BY 8 
ELEMENTS = FHEX2.1 
STEP THETA FR0M 90.0 T0 0.0 BY 4 

NUMBER GPIDS BY 100, ELEMENTS BY 1 
DIVIDE R BY 4 

NUMBER GRIDS BY 10, ELEMENTS BY 10 
PL0T -40., 20. 

ENDM 

$DATA 

TABLE 48,0.,2.,0.,2,4.,0.,4.,2. 


2.3.3 Model: 


Figure 2.3 shows the model generated 
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SAMPLE PROBLEM 2.4 
2.4.1 Description ; 

Generate the solid model corresponding to Problem 1.4. 


2.4.2 MESHLAN Program ; 

M0DEL TEST2.4 
GE0METRY 

SHAPE = FLUID 

B0UNDARY = FUNCTI0N 3 (1.0) 

GRAV =200 
MESH(,1) 

0PR0P =200 
DIVIDE Z BY 3 

NUMBER GRIDS BY 8, ELEMENTS BY 1 
ELEMENTS = FHEX2.1 
DIVIDE THETA BY 8 

NUMBER GRIDS BY 1, ELEMENTS BY 3 
DIVIDE R BY 3 

NUMBER GRIDS BY 100, ELEMENTS BY 100 
PL0T 5., 20. 

ENDM 


2.4.3 Model; 


The three-dimensional plot of the model is shown in Figure 2.4. 
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SAMPLE PROBLEM 2.5 



2.5.1 Description ; 

Model a solid of revolution corresponding to 1.7 using explicitly 
defined axial stations. 


2.5.2 MESHLAN Program ; 

M0DEL TEST2.5 
GE0METRY 

SHAPE - FLUID 
B0UNDARY - TABLE 88 
MESH (.101) 

DIVIDE 2 

NUMBER GRIDS BY 1, ELEMENTS BY 1 
ELEMENTS - FHEX2.1 
STEP THETA FR0M 0.0 T0 90.0 BY 4 

NUMBER GRIDS BY 100, ELEMENTS BY 6 
DIVIDE R BY 3 

NUMBER GRIDS BY 1000, ELEMENTS BY 100 
PL0T -40., 20. 

ENDM 

$DATA 

TABLE 88, 6., 12., 0., 7, 6., 6., 5. 657, 8., 4. 472, 10., 3. 317, 
11., 2., 11. 657,1., 11. 916,0., 12. 


2.5.3 Model; 


Figure 2.5 Illustrates the model generated. 


r 



53 




54 





SAMPLE PROBLEM 3.1 


3.1.1 Description ; 

Generate a three>dlmenslonal model of a 90° segment of a circular 
cylinder r“*5 with a surface tilted 10*. 


3.1.2 MESHLAN Program ; 

M0DEL TEST3.1 
GE0METRY 

SHAPE - GFLUID 

B0UNDARY - FUNCTI0N 4 (3.0) 

GRAV = 100 . 

MESH(,1) 

0PR0P-2OO 
DIVIDE Z BY 3 

NUME-R GRIDS BY 5, ELEMENTS BY 4 
ELEMENTS =• FHEX2,1 
STEP THETA FR0M 0.0 T0 90.0 BY 4 
NUMBER GRIDS BY 1, ELEMENTS BY 1 
DIVIDE R BY 3 

NWffiER GRIDS BY 100, ELEMENTS BY 100 
PL0T -40., 20. 

ENDM 

$DATA 

FUNCTION 4,3.,0.,10.0,0.,0.,0.,1.,0.,5. 


3.1.3 Model; 


Figure 3.1 Illustrates the model generated. 






FIGURE 3.1 
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SAMPLE PROBLEM 3.2 


3.2.1 Description : 

Model a three-dlicenslcnal parabolic solid with a tilted surface, 
a-5*. 


3.2.2 MESHLAN Program : 

M0DEL TEST3.2 
GEOMETRY 

SHAPE - GFLUID 

BOUNDARY - FUNCTION 3 (1.0) 

GRAY - 200 
MESH(,1) 

0PROP - 
DIVIDE Z BY 3 

WjMBER girds by 8, ELEMENTS BY 1 
ELEMENTS - FHEX2,1 
DIVIDE THETA BY 8 

NUMBER GRIDS BY 1, ELEMENTS BY 3 
DIVIDE R BY 3 

NUMBER GRIDS ICO, ELEMENTS BY 100 

PLOT 5., 20. 

ENDM 

$DATA 

FUNCTION 3,10.0,0.0,5.0,5. ,1. ,-.04,0. ,0. ,0. 


3.2.3 Model : 

The t ree-dimenslonal MESHGEN plot Is shown In Figure 3.2. 
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SAMPLE PROBLEM 4.1 


4.1.1 Description ; 

Model the three-dimensional segment of the cylinder defined in 1.1. 
Generate both the tank and fluid models in one step. 


4.1.2 MESHLA:? Program ; 


M0DEL TEST 4.1 
GE0METRY 

SHAPE - TFULL 

B0UNDARY FUKCTI0N 1(3.0) 

GRAY - 20J 
iiESH (1,1001) 

DIVIDE Z BY 3 
SHELL 

NUMBER GRIDS BY 5, ELEMENTS BY 4 
ELEMENTS - QUAD2.1 
PR0PERTY - 100 

S0LID 

NUMBER GRIDS BY 5,ELQIENTS Bf 4 
ELEMENTS “ FHEX2,10i 
PR0PERTY -200 

STEP THETA FR0M 0.0 T0 90.0 BY 4 
SHELL 

NUMBER GRIDS BY 1, ELEMENTS BY 1 
S0LID 

Nm:BER GRIDS BY 1, ELEMENTS BY 1 
DIVIDE R BY 3 

NUMBER GRIDS BY 100, ELEMENTS BY 100 
PL0T -45. ,20. 

LT2 

ENDM 

$DATA 

FUNCTI0N 1,3.,0.,0.,0.,0.,0.,1.,0.,5. 


4.1.3 M odel ; 

Two- and three-dimensional plots for the tank are shown in Figures 
4.1 and 4.i. The three-dimensional fluid is shown in Figure 4.3. 
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SAMPLE PROBLEM 4.2 


4.2.1 Description ; 


Model a section of a tank/fluid system corresponding to the first 
five-inch section of the SRI tank. 


4.2.2 MESHLAN Program ; 

M0DEL TEST4.2 
GE0METRY 

SHAPE - C0NTAINED S0LID 
B0UNDARY - FUNCTI0N 2 (8.66) 

GRAV =200 
MESH (1,1001) 

0PR0P = 100 

STEP Z FR0M 10.0 T0 8.66 BY 2 

NUMBER GRIDS BY 100, ELEMENTS BY 100 
ELEMENTS = QUAD2,1 
STEP Z FR0M 8.66 T0 5.0 BY 4 
SHELL 

NUMBER GRIDS BY 100, ELEMENTS BY 100 
ELEMENTS = QUAD2.201 

f0LID 

NUMBER GRIDS BY 100, ELEMENTS BY 100 
ELEMENTS = FHEX2,100 
PR0PERTY - 100 

DIVIDE R BY 3 

NUMBER GRIDS BY 10, ELEMENTS BY 10 
STEP THETA FR0M 0.0 T0 90.0 BY 4 
SHELL 

NUMBER GRIDS BY 1, ELEMENTS BY 1 
S0LID 

NUMBER GRIDS BY 1, ELEMENTS BY 1 
PL0T -45., 20. 

PLT2 

ENDM 

$DATA 

FUNCTI0N 2,10.,5.,0.,0.,0.,0.,1.,0.,5. 


4.2.3 Model : 


Figures 4.4, 4.5 and 4.6 show the two- and three-dimensional shell 
and three-dimensional fluid, respectively. 
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SAMPLE PROBLEM 5.1 


5.1.1 Description ; 

To Model a 90° sector of a hemisphere defined by a tabular function 
on the Interval z*[1.0t0.0]. The model will be composed of TRIAl 
elements and the shell thickness varies uniformly in the axial 
direction. 


5.1.2 MESHLAN Program ; 

M0DEL TESTS. 0 
GE0METRY 

SHAPE - CAPS 

BOUNDARY = TABLE 36.INTERP0LATE 
MESH(l) 

0PR0P - 100 

STEP Z FROM 1.0 TO 0.0 BY 4 

NUMBER GRIDS BY 100, ELEMENTS BY 100 
ELEMENTS » TRIAl, 1 
THICKNESS VARIES 
STEP THETA FROM 0.0 TO 90.0 BY 4 
NUMBER GRIDS BY 1, ELEMENTS BY 1 
PL0T -45., -20. 

PLT2 

ENDM 

$DATA 

TABLE 36, 1.0, 0.0, 0.0, 9, 0.0, 0.0, 2. 236,. 2, 2. 739,. 3, 3. 162, 
.4, 3. 536,. 5, 3. 873,. 6, 4. 183,. 7, 4. 472,. 8, 5., 1.0 


5.1.3 Model: 


Figures 5.1 and 5.2 show the two- and three-dimensional MESHGEN plots. 
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figure 5.1 
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SAMPLE PROBLEM 5.2 



5.2,1 Deacrlptlon ; 

To model a 120” sector of a hemisphere defined by a tabular function 
on the interval z«[6. 0,12.0]. The model will be composed of TRBSC 
elements with constant properties. 


5.2.2 MESHLAN Program ; 

M0DEL TEST5.2 
GE0METRY 

SHAPE - CAPS 

B0UNDARY - TABLE 88 . INTERP0UTE 
MESH(l) 

0PR0P - 100 
DIVIDE Z BY A 

NUMBER GRIDS BY 1. ELEMENTS BY 100 
ELEMENTS - TRBSC, 1 
STEP THETA FROM 0.0 TO 120.0 BY A 

NUMBER GRIDS BY 100, ELEMENTS BY 1 
PL0T -60., 20. 

PLT2 

$DATA 

TABLE 88, 6., 12., 0. , 7, 6., 6., 5. 657, 8., A. A72, 10., 3. 317, 
11. ,2. ,11.657,1. ,11.716,0. ,12. 


NOTE: When specifying the boundary of a CAPSHELL the Z1 and Z2 values 

on the FUNCTION or TABLE data cards must correspond to the open 
and closed ends of the Idealization, respectively. 


5.2.3 Model ; 

Two- and three-dimensional plots are shown in Figures 5.3 and 5. A. 


69 


i 


1 

1 . 
1 - 
I 


) : 









SAMPLE PROBLEM 5.3 


5.3.1 Description ; 

2 

To model a 90° sector of a shell defined by the parabola 25z«r on 
the Interval [1.0, 0.0]. The model will contain double elements 
(overlapping) of types TRIAI and TRMEM each with different, but 
uniform, properties. The top ring (z»1.0) will be fixed, and sym- 
metric boundaries Imposed at 6-0“ and 6-90“. 


5.3.2 MESHGEN Program ; 

M0DEL TESTS. 3 
GEOMETRY 

SHAPE = CAPS 
B0UNDARY - FUNCTI0N 2 
MESH (101) 

DIVIDE Z BY 6 

NUMBER GRIDS BY 10, ELEMENTS BY 20 
ELEMENTS = TRIAI, 1 
PR0PERTY = 10 
ELEMENTS = TRMEM, 1001 
PR0PERTY = 38 
FIX 123456 AT (1.0) 

STEP THETA FROM 90.0 TO 0.0 BY 6 
NWffiER GRIDS BY 1, ELEMENTS BY 1 
FIX 246 AT (0.0,90.0) 

PL0T -45., -20. 

ENDM 

FUNCTI0N 2,1.0,0.0,0.,0.,1.,-4.,0.,0.,0. 


5.3.3 Model ; 

A three-dimensional structure plot is shown in Figure 5.5. 
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SAMPLE PROBLEM 5.4 


5.4.1 Description ; 

To model a tank closure using the CAPSHELL for the tank defined in 
sample problem 1.8C. 


5.4.2 MESHLAN Program ; 

M0DEL TEST5.4 
GE0METRY 

SHAPE - CAPS 

B0UNDARY = TABLE 73 INTERPOLATE 
MESH (901) 

DIVIDE 2 BY 4 

NU>ffiER GRIDS BY 100, ELEMENTS BY 10 
ELEMENTS = TRIA2.81 
STEP THETA FROM 0.0 TO 90.0 BY 3 
FIX 246 AT (0.0,90.0) 

NUMBER GRIDS BY 10, ELEMENTS BY 1 
PL0T - 45 ., - 20 . 

ENDM 

$DATA 

TABLE 73,3. ,0. ,0., 7,0., 0.,1. ,.3, 2.,. 7, 3. ,1.1,4. ,1.6, 
5. ,2.25,6. ,3. 


5.4.3 Model; 


Figure 5.6 shows a three-dimensional MESHGEN plot of this model 





SAMPLE PROBLEM 6.1 


6.1.1 Description ; 

2 2 

Model a quadrant of the hemisphere defined by z +r = 25 on the 
interval [.,-5.]. Use H'JEDGE elements. 


6.1.2 MESHLAN Program : 

M0DEL TEST6.1 
GE0METRY 

SHAPE =■ CAPFLUID 
B0UNDARY = FUNCTI0N 8 
MESH(,1) 

DIVIDE Z BY 4 

NUMBER GRIDS BY 100, ELEMENTS BY 100 
ELEMENTS = FWEDGE,! 

STEP THETA FR0M 0.0 T0 90.0 BY 2 
NU>IBER GRIDS BY 1, ELEMENTS BY 1 
DIVIDE R 

NUMBER GRIDS BY 10, ELEMENTS BY 10 
PL0T -4O.,2O.,N0NUM 
ENDM 
$DATA 

FUNCTI0N 8,0.,-5.,0.,l.,0.,l.,0.,5.,0. 


NOTE: The command DIVIDE R does not need an increment since the 

number of R steps must equal the z division. The value of 
zl in the FUNCTION (or TABLE) definition must be the open 
end of the solid. CAPFLUID SHAPES may not produce plots via 
the PLT2 option. 


6.1.3 Model: 


Figure 6.1 illustrates the three-dimensional model. 





SAMPLE PROBLEM 6.2 
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6.2.1 Description ; 

2 

Model half of the body generated by the parabola 25z = r on the 
interval z>«[4.0,0.]. Use FHEX2 elements. 


6.2.2 MESHLAN Program ; 

M0DEL TEST6.2 
GE0METRY 

SHAPE - CAPF 
B0UNDARY « FUNCTI0N 1 
MESH(,1) 

DIVIDE Z BY 4 

NUMBER GRIDS BY 50, ELEMENTS BY 50 
ELEMENTS - FHEX2.1 
STEP THETA FR0M 90.0 T0 0.0 BY 2 
NUMBER GRIDS BY 1, ELEMENTS BY 1 
DIVIDE R 

NUMBER GRIDS BY 10, ELEMENTS BY 10 
PL0T -40.0,-20.0 
ENDM 
$DATA 

FUNCTI0N 1,4. 0,0. 0,0. ,0. ,25. ,-l. ,0. ,0. ,0. 


6.2.3 Model; 


Figure 6.2 shows the three-dimensional model. 
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SAMPLE PROBLEM 6.3 
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6.3.1 Description ; 

Generate a model for a hemisphere defined by z^ + r^-25 cn the 
Interval [0.0, 5.0]. Use FWEDGE elements. 


6.3.2 MESHLAN Program ; 

M0DEL TEST6.3 
GE0METRY 

SHAPE - CAPF 
B0UNDARY - FUNCTI0N 3 
MESH(,1) 

DIVIDE Z BY 5 

NUMBER GRIDS BY 100, ELEMENTS BY 100 
ELEMENTS - FWEDGE,! 

DIVIDE THETA BV 8 

NUMBER GRIDS BY 1, ELEMENTS BY 1 
DIVIDE R 

NUMBER GRIDS BY 10, ELEMENTS BY 10 
PL0T -5.,-2O.,N0NUM 
PLT2 
ENDM 
$DATA 

FUNCTI0N 3,0.,5.,0.,1.,0.,1.,0.,5.,0. 


6.3.3 Model; 


A three-dimensional MESHGEN plot is shown in Figure 6.3. 
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SAMPLE PROBLM 7.1 
7.1.1 Description ; 

2 2 

Generate a model for the upper hemisphere of the sphere z -f r "25 
on z»[0.t5.]. Include both fluid and structure models. 


7.1.2 MESHLAN Program ; 

M0DEL TEST7.1 
GE0MEIRY 

SHAPE - CAPB 
B0UNDARY - FUNCTI0N 1 
GRAVITY - 1 
MESH (1,1001) 

DIVIDE Z BY 4 
SHELL 

NUMBER GRIDS BY 100, ELEMENTS BY 100 
ELEMENTS «> TRMEM,1 
PR0PERTY - 100 

S0LID 

NUMBER GRIDS BY 100, ELEMENTS BY 100 
ELEMENTS - FWEDGE,1001 
PR0PERTY - 200 

STEP TEHTA FROM 0.0 TO 90.0 BY 4 
SHELL 

NUMBER GRIDS BY 1, ELEMENTS BY 1 
S0LID 

NUMBER GRIDS BY 1, ELEMENTS BY 1 
DIVIDE R 

NUMBER GRIDS BY 10, ELEMENTS BY 10 
PL0T -45., 20. 

PLT2 

ENDM 

$DATA 

FUNCTION 1,0.,5.,0.,1.,0.,1.,0.,5.,0. 


7.1.3 Model ; 

The two- and three-dimensional MESHGEN plots of this model are shown 
In Figures 7.1 - 7.3. 
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SAMPLE PROBLEM 7.2 


7.2.1 Description ; 

Generate a model for the same function as problem 7.1 on the 
interval [0,-5.]. 


7.2.2 MESHLAN Program ; 

M0DEL TEST7.2 
GE0METRY 

SHAPE » CAPB 
BOUNDARY = FUNCTI0N 2 
GRAVITY - 225 
MESH (1,1001) 

DIVIDE Z BY 4 
SHELL 

NUMBER GRIDS BY 100, ELEMENTS BY 100 
ELEMENTS = TRIA1,1 
PR0PERTY = 100 

S0LID 

NUMBER GRIDS BY 100, ELEMENTS BY 100 
ELEMENTS = FWEDGE,1001 
PR0PERTY =200 

STEP THETA FROM 0.0 TO 90.0 BY 4 
SHELL 

NUMBER GRIDS BY 1, ELEMENTS BY 1 
S0LID 

NUMBER GRIDS BY 1, ELEMENTS BY 1 

DIVIDE $ 

NUMBER GRIDS BY 10, ELEMENTS BY 10 
PL0T -45., 20. 

PLT2 

ENDM 

$DATA 

FUNCTI0N 2,0.,-5.,0.,l.,0.,l.,0.,5.,0. 


NOTE; For a CAPB SHAPE option, the value z, on the FUNCTION card 
(or TABLE) must define the open end of the figure 


7.2.3 Model ; 

Figures 7.4 - 7.5 .‘jhow the two- and three-dimensional MESHGEN plots. 
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SAMPLE PROBLEM 8.1 


8.1.1 Description ! 

To mode? a 15° segment of the SRI tank. The tank Is composed of 
two sections, the top 5 inches ai.! a cylinder, r=>5, and the last 
5 Inches are a hemisphere. The upper portion of the tank vill be 
modf^led using the TFULL shape. The bottom section will use TANK 
and CAPFLUID shapes. Figure 8.1 Illustrates the model desired. 


8.1.2 MESHLAN Program : 

H0DEL SRITANKl 
GE0METRY 

SHAPE = TFULL 

B0UNDARY = FUNCTI0M 1 (8.66) 

GRAY = 200 
MESH (5010, 1012) 

STEP Z FRtf'M 10.0 T0 8.66 BY 2 

NUMBER GRIDS BY 1, ELEMENTS BY 1 
FIX 123456 AT (10.0) 

ELEMENTS = QUAD2,1 
PR0PERTY » 100 
STEP Z FR0M 8.66 T0 5.0 BY 4 
SHELL 

NUMBER GRIDS BY 1, ELEMENTS BY 1 
ELEMENTS » QUAD2,3 
PR0PERTY = 100 

S0LID 

NUMBER GRIDS BY 1, ELEMENTS BY 100 
ELEMENTS = FHEX2,301 
PR0PERTY = 200 

STEP THETA FR0M 0.0 T0 15.0 BY 1 
SHELL 

NUMBER GRIDS BY 1000, ELEMENTS BY 1 
S0LID 

NUMBER GRIDS BY 1000, ELEMENTS BY 1 
DIVIDE R BY 8 

NUMBER GRIDS BY 1-0, ELEMENTS BY 1 
PL0T 180., 10. 

ENDM 

M0DEL SRITANK2 
GE0METRY 

SHAPE = TFULL 
B0UNDARY = FUNCTI0N 2 
GRAVITY =200 
MESH(5016,1016) 

DIVIDE Z BY 5 


SHELL 

NUMBER GRIDS BY 1, ELEMENTS BY 1 
ELEMENTS - QUAD2.7 
PROPERTY - 100 

S0LID 

NUMBER GRIDS BY 1, ELEMENTS BY 100 
ELEMENTS - FHEX2.701 
PR0PERTY - 200 

STEP THETA FR0M 0.0 T0 15.0 BY 1 
SHELL 

NUMBER GRIDS BY 1000, ELEMENTS BY 1 
S0LID 

NUMBER GRIDS BY 1000, ELEMENTS BY 1 
DIVIDE R BY 8 

NUMBER GRIDS BY 100, ELEMENTS BY 1 
PL0T 180. ,10. 

ENDM 

M0DEL SRITANK3, SAVE, NEW 
6E0METRY 

SHAPE - CAPF 
B0UNDARY =• FUNCTI0N 3 
MESH(,1021) 

DIVIDE Z BY 8 

NUMBER GRIDS BY 1, ELEMENTS BY 100 
ELEMENTS = FHEX2,1201 
PR0PERTY = 200 

STEP THETA FR(»I 0.0 TO 15. BY 1 

NUMBER GRIDS BY 1000,ELEMENTS BY 1 
DIVIDE R 

NUMBER GRIDS BY 100, ELEMENTS BY 1 
PL0T 180. ,-10. 

ENDM 

M0DEL SRITMK4,SAVE 
GE0METRY 

SHAPE =" TANK 
B0UNDARY » TABLE 50 
MESH (5021) 

DIVIDE Z 

NUMBER GRIDS BY 1, ELEMENTS BY 1 
ELEMENTS = QUAD2,12 
PR0PERTY - 100 

STEP THETA FR0M 0.0 T0 15.0 BY 1 

NUMBER GRIDS BY 1000, ELEMENTS BY 1 
PL0T 180., 10. 

ENDM 

M0DEL SPRITANK5 
GRAVITY = 200 

FIND B0UNDARIES SRITANK4 , SRITANK3 

PUNCH 

ENDM 


NOTES; 


a. SRITANK3 is to be saved on a SAVE file declared 
NEW In order to generate CFLSTR Bulk Data later. 

b. SRITANK4 Is added to the SAVE file. 

c. SRITANK5 retrieves the structure (SRITANK4) and 
fluid (SRITANK3) models from the SAVE file and 
generates CFLSTR cards for their union. 


8.1.3 Model ; 

Figure 8.1 defines the model desired. The MESHGEN plots generated, 
as defined in the MESHLAN programs above, are shown in Figures 
8 . 2 - 8 . 7 . 
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